Topological Classification of Gapped Spin Chains :Quantized Berry Phase as a Local 
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We characterize several phases of gapped spin systems by local order parameters defined by 
quantized Berry phases^. This characterization is topologically stable against any small perturbation 
as long as the energy gap remains finite. The models we pick up are S = 1, 2 dimerized Heisenberg 
chains and S = 2 Heisenberg chains with uniaxial single-ion-type anisotropy. Analytically we also 
evaluate the topological local order parameters for the generalized Affieck-Kennedy-Lieb-Tasaki 
(AKLT) model. The relation between the present Berry phases and the fractionalization in the 
integer spin chains are discussed as well. 



PACS numbers: 75.10.Jm, 03.65.Vf, 73.43.Nq, 75.10.Pq 



I. INTRODUCTION 

Characterizing quantum many-body systems is one of 
the important topics in condensed matter physics. The 
Ginzburg-Landau (GL) theory has been quite successful 
to describe many phases based on a concept of the sym- 
metry breaking and the local order parameter. Despite 
its remarkable success, novel types of phases which are 
not well described by the (classical) local order parame- 
ters have been found in many systems. Concepts of topo- 
logical order and quantum order are trial to overcome 
the difficulties of the classical GL theory with the sym- 
metry breaking^''^''*. One of the characteristic features 
of the topological insulators is that localized states, such 
as the edge states, appear near the system boundaries, 
even though the system without boundaries has a finite 
energy gap. Examples of such systems are quantum Hall 
liquida^i^ii, Haldane spin systema^i^, polyacetyleneiS, 
and spin-Peierls systeroii. Recently it has become clear 
that the bulk-edge correspondence^ii has an intimate 
relation to the entanglement entrop y^^i^^'^"^ , which has 
been discussed to detect non trivial structures of topo- 
logically ordered stateai^ii^ii^ii^. 

Recently one of the authors proposed to use another 
quantum quantity, quantized Berry phases^i, to define 
a topological local order parameteri^iii. One can define 
a topological local order parameter by the Berry phases 
even though there is no classical order parameter. The 
Berry phases is a typical quantum quantity based on the 
Berry connection which is defined by the overlap between 
the two states with infinitesimal difference. It implies 
that the topological local order parameter defined ini is 
a quantum order parameter that does not have any cor- 
responding classical analogs. Further it has a conceptual 
advantage for the topologically ordered phases, since it 
is quantized to or tt (mod 27r) when the ground state 
is invariant under some anti-unitary transformation. It 
implies a topological stability that the quantized Berry 
phase does not change against any small perturbation. 



The Berry phase is given by an integration of the Berry 
connection defined by the local U{1) twist on a link of 
a lattice. Then the quantum phases can be categorized 
by the texture pattern of the Berry phases (0 or tt). It 
has been successfully applied to several gapped quantum 
systems. For example, the ground states of the S — 1/2 
dimerized Heisenberg models (in one and two dimensions 
even with frustrations )i can be characterized by the pat- 
tern of TT Berry phases on the bonds which indicate the 
locations of dimer singlets. In a case of the t-J modeli^, it 
is characterized by the texture pattern of the non-Abelian 
Berry phase, which describes itinerant singlets. Also, for 
the S* = 1 Heisenberg model, its ground state as the Hal- 
dane phase was characterized by the uniform tt Berry 
phases. This topological order parameter also clearly de- 
scribe a quantum phase transition between the Haldane 
phase and the large-D phasesi^. 



In this paper, we calculate the topological local or- 
der parameter by the quantized Berry phase for several 
gapped quantum spin chains. There are substantial num- 
bers of studies for the Haldane phas o^°i^^'^^i^^ . Then it 
has been clarified that the Haldane phase can be char- 
acterized by the hidden Z2 x Z2 symmetry breaking^iiS^ 
which describes a non-locality of the Haldane phase (by 
the string order parameters On the other hand, 
the topological order parameter by the Berry phases is 
local and quite useful for the S — 1 case to describe 
the phase and the quantum phase transitioni^-. Here we 
further investigate generic situations, such as the several 
Haldane phases in the dimerized S = 1,2 Heisenberg 
chainsSSi^ and the S = 2 Heisenberg chain with uniaxial 
single-ion-type anisotrop y^^i'^^i^'^i'^^ . We also study the 
Berry phase of the generalized valence-bond-solid (VBS) 
state analytically and interpret the numerical results in 
terms of the reconstruction of the valence-bonds. 



II. DEFINITION OF THE BERRY PHASE 



Let us start with defining the Berry phase in a quan- 
tum spin system. The Berry phase is defined when the 
Hamiltonian has parameters with periodicity assuming a 
finite energy gap between the ground state and the ex- 
cited statesii. For the parameter dependent Hamiltonian 
H{(f)), the Berry phase 7 of the ground state is defined 
as 



h = f A{<j))d<j), (1) 

^0 



where A{(f) is the Abelian Berry connection obtained 
by the single- valued normalized ground state |GS (</>)) of 
H{4>) as A{4)) = (GS(0)|a0|GS(0)). This Berry phase is 
real and quantized to or tt (mod 27r) if the Hamiltonian 
H{(j)) is invariant under the anti-unitary operation 0, i.e. 
[H{<j)), e] = 0^. Note that the Berry phase is "undefined" 
if the gap between the ground state and the excited states 
vanishes while varying the parameter (j). We use a local 
spin twist on a link as a generic parameter in the defini- 
tion of the Berry phase^. Under this local spin twist, the 
following term S~ + in the Hamiltonian is re- 

placed with e"f'S+Sj+e-"l'S-S+, where 5*=^ = Sf±iSf. 
The Berry phase defined by the response to the local spin 
twists extracts a local structure of the quantum system. 
By this quantized Berry phase, one can define a link- 
variable. Then each link has one of the three labels: "0 
bond" , "tt bond" , or "undefined" . It has a remarkable 
property that the Berry phase has topological robustness 
against the small perturbations unless the energy gap be- 
tween the ground state and the excited states closes. On 
the other hand, the "undefined" indicates an existence of 
the quantum phase transition. In order to calculate the 
Berry phase numerically, we introduce a gauge-invariant 
Berry phasoii^ on a lattice. It is defined by discretizing 
the parameter space of into N points as 



In = - ^ argAAr(0„), 0„ = — n, (2) 

n— 1 



where ((/)„) is defined by Ais[{(j)n) = 

(GS((/)„)|GS((/)„+i)), (j)N+i = (jJi- We expect 
7 = limjv_,c» 7Ar. To calculate -fN, we use the 
Lanczos method to diagonalize the Hamiltonian in the 
subspace oi^^Sf =0. 
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FIG. 1: (Color Online) The Berry phases on tlie local link of 
(a) the S = 1 periodic N = 14 and (b) the 5 = 2 periodic 

= 10 dimerized Heisenberg chains. The Berry phase is 
TT on the bold line while that is on the other line. The 
phase boundaries in the finite size system are 9ci = 0.531237, 

= 0.287453 and 0^3 = 0.609305, respectively. The Berry 
phase in (a) and (b) has an inversion symmetry with respect 
to S = 7r/4. A schematic VBS picture of the ground state 
is assigned to each phase. Dots, bold lines, and open circles 
denote the S = 1/2, singlet dimers, and the operations of 
symmetrization, respectively. 



III. 5 = 1,2 DIMERIZED HEISENBERG 
MODELS AND 5 = 2 HEISENBERG MODEL 
WITH UNIAXIAL ISOTROPY 

A. 5 = 1,2 dimerized Heisenberg models 

First we consider S — 1, 2 dimerized Heisenberg mod- 
els 

N/2 

H — {JlS2i ■ 821+1 + J2S2i+l ■ 821+2) , (3) 

1=1 

where Si is the spin-1 or 2 operators on the i-th site and 
N is the total number of sites. The periodic boundary 
condition is imposed as S^+i = Si for all of the models 
in this paper. Ji and J2 are parameterized as ,/i = sm9 
and J2 — cos9, respectively. We consider the case of 
< 6 < tt/2 in this paper. The ground state is composed 
of an ensemble of A''/2 singlet pairs in limits of 6* ^ and 
9 — > 7r/2. The system is equivalent to the isotropic anti- 
ferromagnetic Heisenberg chain at 6* = 7r/4. Based on the 
VBS picture, we expect a reconstruction of the valence 
bonds by changing 9. 

Figures. [BJa) and (b) show the 9 dependence of the 
Berry phase on the link with Ji coupling and J2 cou- 
pling with S' = 1, TV = 14 and S* = 2, TV = 10, respec- 
tively. The region with the Berry phase tt is shown by 
the bold line. There are several quantum phase transi- 
tions characterized by the Berry phase as the topologi- 
cal order parameters. The boundary of the two regions 
with different Berry phases and tt does not have a well- 
defined Berry phase, since the energy gap closes during 
the change of the local twist parameter cj). Since the 
Berry phase is undefined at the boundaries, there exists 
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FIG. 2: (Color Online) The Berry phases on the local link of 
the S — 2 periodic A'' = 10 Heisenberg chain with single-ion 
anisotropy. The notations are the same as the Fig[T] An up 
(down) triangle denotes an up (down) spin-1/2. 



the level crossing which imphes the existence of the gap- 
less excitation in the thermodynamic limit. This result is 
consistent with the results previously discussed'^'', which 
the general integer-S* extended string order parameters 
changes as the dimerization changes. The phase diagram 
defined by our topological order parameter is consistent 
with the one by the non-local string order parameter. In 
an = 10 system with S = 2, the phase boundaries 
are 6c2 — 0.287453, Ocs — 0.609305, and it is consistent 
with the results obtained by using the level spectroscopy 
which is based on conformal field theory techniques^. 
Especially in the one dimensional case, the energy di- 
agram of the system with twisted link is proportional 
to that of the system with twisted boundary conditions. 
However, our analysis focuses on the quantum property 
of the wave functions rather than the energy diagram. 



B. S = 2 Heisenberg model with uniaxial anisotropy 

As for the = 2 Heisenberg model with Z3-term, we 
use the Hamiltonian 



N 

H = Y.[jS^,-S,+,+D{StY 



(4) 



Figure. [5] shows the Berry phase of the local link in the 
S = 2 Heisenberg model + Z?-term with N=10. The 
parameter J = 1 in our calculations. The region of the 
bold line has the Berry phase tt and the other region 
has the vanishing Berry phase. This result also makes 
us possible to consider the Berry phase as a local order 
parameter of the Haldane spin chains. Our numerical 
results for finite size systems support the presence of the 
intermediate _D-phase^. 



IV. INTERPRETATION OF THE NUMERICAL 
RESULTS BY VALENCE-BOND SOLID STATE 

Let us now interpret our numerical results in terms 
of the VBS state picture. The VBS state is the exact 
ground state of the Affleck-Kennedy-Lieb-Tasaki(AKLT) 
model^. We shall calculate the Berry phase of the gener- 
alized VBS state with the aid of the chiral AKLT model^^ 
and its exact ground state wave function. The chiral 
AKLT model is obtained by applying 0(2) rotation of 



spin operators in the original AKLT model. In our cal- 
culation, it is convenient to introduce the Schwinger bo- 
son representation of the spin operators as = ajbi, 



Qib], and Sf 



blbi)/2. ai and hi sat- 
isfy the commutation relation [ai,aj] = [^i,^]] = Sij 
with all other commutators vanishing"^^. The constraint 
a|ai + blbi — 2Si is imposed to reproduce the dimen- 
sion of the spin Si Hilbert space at each site. In general, 
the ground state of the chiral AKLT model having Bij 
valence bonds on the link {ij) is written as 



'' J 



I vac) (5) 



38 



This state has nonzero average of vector spin chirality 



{Si X Sj ■ z) unless the twist parameter 



or TT. 



This state is a zero-energy ground state of the following 
Hamiltonian: 



N J„ 



(6) 



where Jmax = + 2_Bi^i+i + i3i+i^i+2) /2, Jmin = 

{Bi-i^i -\- -Bi_|_i j_|_2) /2-1-1 and Aj is the arbitrary positive 
coefficient. iVi-t-iP] is the polynomial in Si ■ Si+i and 
act as a projection operator projecting the bond spin 
Ji.i+i = Si + Si+i onto the subspace of spin magnitude 
J. The replacement 



i+l c- c + 

'-'i '-'i+l 



(7) 

in Si ■ Si+i produces Pf^^^j^[(j)i^i+i] in Eq. (0). 

Now we shall explicitly show that the Berry phase of 
the VBS state extracts the local number of the valence 
bonds Bij as i?.y7r(mod 2tt). Let us now consider the 
local twist of the parameters 4'ij — (t>Sij^i2 and rewrite 
the ground state as To calculate the Berry 

phase of the VBS state, the following relation is useful: 



«7i2 = iBi2Tr + i 



Im[(0|9^|0)]/AA(0)d<^, (8) 



where 712 is the Berry phase of the bond (12) and 
^/{(j)) = {(t>\4>)- Note that the first term of the right hand 
side comes from the gauge fixing of the multi- valued wave 
function to the single-valued function. Then, the only 
thing to do is to evaluate the imaginary part of the con- 
nection 

Let us first consider the 5* = 1 VBS state as the sim- 
plest example. In this case, -B^.i+i = 1 for any bond and 
the VBS state with a local twist is given by 

N 

(9) 

We impose the periodic boundary condition, i.e., Oat+i = 
fli and 6Ar+i = 61. It is convenient to introduce the 
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singlet creation operator s''' = {a\bl—b\al) and the triplet 
{Jz = 0) creation operator = {a\bl + b\al). We can 
rewrite the bond (12) part of the VBS state (e'-^/^^t^t _ 
e-^'^^%\al) as (cos|st + jsin|tt). Then |(/)) and 9^|0) 
can be written as 



A) =cos-|0) + isin-|l), 

1 </> I , i 4> I \ 
--sm-|0) + -cos-|l), 



(10) 



where 



AT 



|0) = s^Ylialbl^-blal,)\vac), 

i=2 
N 

|1) = ^^n(«I^>i-^H+i)|vac). 



It is now obvious that the imaginary part of {4>\d^\4>) 
vanishes since the state |1) having a total spin S'totai = 1 
is orthogonal to the state |0) with S'totai = 0. Therefore, 
the Berry phase of this state is given by 712 = tt. Next 
we shall consider a more general situation with arbitrary 
Bij. We can also express the VBS state with a local twist 
on the bond (12) in terms of and as 



cos— s' +«sm— r' 
^22 



B12 



N 



Y[ialbl,-blal,f'-^^\y^c) 



i=2 



By using the binomial expansion, can be rewritten as 



k=0 



i:(--')(c»|)'-(M„|)|.,, 



where |fc) = (st)Si2-fc(tt)fe(. . . )|vac) is the state with k 
triplet bonds on the link (12). (• • ■ ) denotes the rest of 
the VBS state. In a parallel way. 



k=0 ^ ^ 

X (^isin^)''(fccot^ - (S12 -fc)tan^)|fc). 

To see that the imaginary part of {(\)\d^\(\)) is zero, we 
have to show that Im(fc|Z) — when k and I have the 
same parity(even or odd) and Re(fc|Z) = when k and 
I have different parities. This can be easily shown by 
using the coherent state representation of the Schwinger 
bosons (see APPENDIX). Then using the relation (H]), 
we can obtain the Berry phase as 

7u = -SyTT, (mod 27r). 

This result means that the Berry phase of the general- 
ized VBS state counts the number of the valence-bonds 
on the bond (ij). One valence-bond has the tt Berry 



Bi2-k 



phase. Finally, it should be stressed that our calculation 
of the Berry phase is not restricted to one-dimensional 
VBS states but can be generalized to the VBS state on 
an arbitrary graph^^ as long as there is a gap while vary- 
ing the twist parameter. 

Now, let us consider the previous two models in terms 
of the VBS picture. For the 5 = 2 dimerized Heisen- 
berg model, the number of the valence bonds changes as 
the changes (see FigllJ. Since the number of the va- 
lence bonds on a local link can be computed by the Berry 
phase, we can clearly see that the reconstruction of the 
valence bonds occurs during the change of the dimeriza- 
tion. Thus, the result of the Berry phase is consistent 
with the VBS picture. For the S — 2 Heisenberg chain 
with single-ion anisotropy, the valence bonds are broken 
one by one as D increases as we can see in the Fig. [5] 
We see that the Berry phase reflects the number of the 
valence bonds as well as the previous dimerized Heisen- 
berg chain. This can be understood as a fractionalization 
since the basic objects of the present integer spin chains 
are spin- 1/2 singlets. 



V. THE RELATION BETWEEN THE BERRY 
PHASE AND THE ENTANGLEMENT ENTROPY 



Moreover, the Berry phase of generahzed VBS state re- 
lates to the number of the edge states which emerge when 
the spin chain has edges-. Thus, it detects the property 
of the topological phase. Since the entanglement entropy 
also detects such phasesi^ii^i^, we clarify the relation 
between the Berry phase and the entanglement entropy. 
The entanglement entropy of our generalized VBS state 
in thermodynamic limit is Sa ~ J2{ij)£dA^'^si^v + 
where dA denotes the set of the bonds on the boundary 
of subsystem A. It counts the number of the edge states 
^odgc as Sa = log^cdgo- Thus, the Berry phase is related 
to the entanglement entropy in generalized VBS states 
via the edge states in the thermodynamic limit^. 



VI. CONCLUSION 

In conclusion, we have shown that the topological lo- 
cal order parameter defined by quantized Berry phases is 
useful to classify the phases of various spin chains such as 
the Haldane phase. In our calculations, the Berry phase 
is locally defined and does not need nonlocal calculations. 
It is also useful to estimate the order parameter from the 
finite size systems since it is quantized even in the fi- 
nite size systems. The property of the phase is revealed 
in terms of the texture pattern of the Berry phase. We 
have also analytically studied the Berry phase of the gen- 
eralized VBS state and found that the Berry phase picks 
up the number of singlets on the local link. 
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APPENDIX 

In this appendix we show that Re(A:|^) = when k and 
I have different parities and Im(fc|/) ~ when k and I 
have the same parity by simple symmetry arguments. To 
show them, it is convenient to introduce a spin coherent 
stated. For a point Cl = (sinScost/), sin6'sin0, cos9) on the 
unit sphere, the spin coherent state at each site is defined 
as 



= , — vac) 



(13) 



where {u,v) = {cos{e/2)e"l'/'^,sm{0/2)e~''*'/^) are spinor 
coordinates. Using jfJ), the resolution of the identity is 
given by 



OC I 1 

/ = —— I dn\n}(n\, 

47r 



(14) 



i^(f^l,f^2) 



-(cos^i 



cos6'2 — isin6'isin6'2sin((/)i — </)2))- 



Here we have already used the following relation: 



(vac a 



s-i 



b^+'\n) = ^/(2sy.u 



S-l^S+l^ 



First we consider 
the case where k and / have different parities. In this 
case, k — l is odd and hence K{rii,ri2)''~^ changes its sign 
under the change of variables {9j,4>j) to (tt — Oj,—(pj) 
{j = 1,2, N). Since the other part of the integrand 
is invariant under this change of variables, we obtain 
{k\l) ~ 0. Therefore it is now obvious that Re(fc|/} ~ 
when k and I have different parities. Next we consider 
the case where k and I have the same parity. In this case, 
k — l is even. Thus we set k — l~ 2m (m G N) and expand 

K{ni,n2f"' as 



m,,^i2?^ = (i) "^ r™ (COS01-COS02)''"-" 

ri=0 ^ ^ 

X (-i)"(sin6iisin6l2sin((/)i -(/)2))". (16) 



where / denotes a (25' + l)-dimensional identity matrix. 
Let us now consider the inner product {k\l). We can set 
k > I without loss of generality. Inserting the resolu- 
tion of the identity (fT4|) between (fc| and |^), the integral 
representation of the inner product can be obtained as 



COS01COS02 ) K{Qi,Q2) 



k-l 



N 

n 



1 - fii • 



(15) 



The imaginary part of K{fli, 172)^™ comes from the con- 
tribution of the odd n's in the above summation. Now 
we consider the following change of variables: {9j,(f>j) 
to {6j,—(j)j) {j = 1,2,...,N). Under this transforma- 
tion, lm[K{Cli,Ct2)'^"^] changes its sign. On the other 
hand, the other part of the integrand in Eq. (|15p is real 
and invariant under this change of variables. Therefore, 
Im(fc|Z) = when k and I have the same parity. Finally, 
we remark that the generalization of the above result to 
the VBS state on an arbitrary graph is almost trivial 
since we have not used a specific property of the one- 
dimensional VBS state in our proof. 
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